Density perturbations in braneworld cosmology and primordial black holes 
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We study, by numerical methods, the time evolution of scalar perturbations in radiation era of 
Randall-Sundrum braneworld cosmology. Our results confirm an existence of the enhancement of 
perturbation amplitudes (near horizon crossing), discovered recently. We suggest the approximate 
solution of equations of the perturbation theory in the high energy regime, which predicts that the 
enhancement factor is asymptotically constant, as a function of scale. We discuss the application of 
this result for the problem of primordial black hole production in braneworld cosmology. 
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I. INTRODUCTION 



During last decade, braneworld cosmological scenar- 
ios, in which our 4D Universe is realized as a hypersur- 
face embedded in a higher-dimensional spacetime have 
attracted much attention. In first scenarios of this kind, 
suggested as early as in 1980's [H, Hj, it had been shown 
that matter fields can be confined to a field-theoretical 
domain wall (topological defect) in a world with non- 
compact extra dimensions. The progress in string theory 
in subsequent years, especially the discovery of D-branes, 
has revived interest to the idea of braneworlds. In gen- 
eral, the string theory is quite promising, it may provide 
an unified description of gauge interactions and gravity. 
In the present context, it is most important that it pre- 
dicts the existence of p-branes, {p + l)-dimensional sub- 
manifolds of the 10 (or 11) - dimensional spacetime on 
which open strings end. Gauge particles and fermions 
which correspond to string end points can only move 
along these p-branes, while gravitons can propagate in 
the full spacetime ("bulk"). It is tempting to assume 
that our (3 + l)-dimensional spacetime is such a 3-brane. 
If only gravity can probe the bulk, the extra dimensions 
can be very large (in comparison with the smallest length 
scale tested, so far, in particle physics, ^ 10^^^ cm). It 
had been assumed in [3|, that the extra dimensions are 
compact, in analogy with the old Kaluza-Klein (KK) pic- 
ture Slightly later, in works by Randall and Sundrum 
@i @| , it was pointed out that this condition is not neces- 
sary and the extra dimension may be even non-compact. 

The Randall-Sundrum (RS) model is of particular in- 
terest due to its relative simplicity, in spite of the fact 
that it includes nontrivial gravitational dynamics. In the 
RS2 model [1| a single brane is embedded in a anti - 
de Sitter (AdS) bulk and, although the 5th dimension 
extends infinitely, the warped structure of the bulk ge- 
ometry (i.e., the curvature of the bulk spacetime) leads 
to a recovery of the standard General Relativity (GR) on 
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the brane at scales larger than the bulk curvature scale 
£. In particular, Newton's law is recovered at large dis- 
tances and the Friedmann's equation for the evolution of 
the Universe is obtained at low energy. 

At high energies, i.e., in the very early Universe, the 
Friedmann equation differs substantially from GR by a 
correction term which is proportional to p/cr, where p is 
the density of brane matter, and a is the brane tension. 
This term leads to a faster Hubble expansion at high 
energies. Inflationary expansion of the Universe is also 
modified in brane cosmology: the evolution of the infla- 
ton field is more strongly damped, and the brane Uni- 
verse inflates at much faster rate than what is expected 
from standard cosmology. Another important effect at 
high energies is the excitation of KK-modes which escape 
from our brane into the 5D bulk, leading, in particular, 
to the suppression of the power spectrum of inflationary 
gravitational wave background. 

Cosmological perturbation theory in braneworld cos- 
mology also has some distinct features The equa- 
tions of the perturbation theory contain high-energy cor- 
rections (^ p/c) similar to those in the Friedmann equa- 
tion and, in addition, the correction terms arising from 
the fluctuations of the bulk geometry. Perturbations on 
brane, e.g., the scalar perturbations (which we are inter- 
ested in) are coupled with the bulk perturbations. Tech- 
nically, in a case of the scalar perturbations and AdS 
bulk, the problem is reduced to the solution of a system 
of equations for the density contrast variable and the 
so-called master variable (it appears that all quantities 
describing the bulkjperturbations are written in terms of 
this variable 0,[l[l3i]). 

In the context of braneworld models, a question about 
existence and evolution laws of the higher-dimensional 
black holes is very interesting and important. In a model 
with the 5th large extra dimension, a physically meaning- 
ful black hole solution is the 5D-Schwarzschild [13, [l^ ; 
if the horizon size is sufficiently small compared with an 
effective size of the extra dimension. Really, it is natural 
to assume that primordial braneworld black holes formed 
in the early Universe with a horizon size <C ^ would be 
described by a 5D Schwarzschild metric because in this 
case the AdS curvature has very little effect on the ge- 
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ometry. Numerical calculations support the existence of 
static solutions for such small [1^. However, the re- 
sults of these calculations cannot be extrapolated to the 
case Vs ~ £. 

Unfortunately, an exact solution representing a local- 
ized and stable black hole is known only in 4D braneworld 
model , whereas the corresponding solution in the 5D 
braneworld model has not been found. The process of 
the gravitational collapse on the brane is very compli- 
cate, due to, in particular, gravitational interaction be- 
tween the brane and the bulk (see, e.g., [l§]). Even in the 
simplest case of RS- type brane, and Oppenheimer-Snyder 
(OS) - like collapse, braneworld gravity introduces impor- 
tant new features in the black hole formation process (the 
high energy- and KK-corrections to the field equations of 
GR, i.e., the same corrections which affect the expansion 
of the early Universe, are also efficient here). These fea- 
tures lead to a non-static exterior of the black hole [l^ 
in the case of the OS-collapse. Moreover, there are argu- 
ments [13:1m based on AdS/CFT-correspondence, that 
the non-static behavior exists also in a general collapse. 
If, really, the black hole solutions in braneworld scenarios, 
for a black hole larger than AdS radius, are quite different 
from those in 4D GR (i.e., if, as authors of [l^, [2l| ar- 
gue, these solutions are necessarily non-static and predict 
short lifetime of large black holes due to the strongly en- 
hanced evaporation) , there is an unique chance to probe 
the extra dimension by astronomical observations of mas- 
sive black holes. 

Predictions for an evolution of the small (r^ <C i) black 
holes are less dramatic (and less speculative) . The differ- 
ences from the 4D case are reduced to a larger probability 
of accretion, in the high energy regime (due to the fact 
that in this regime the radiation density is proportional 
to t^^ rather than i~^) and to a relative increase of the 
primordial black hole (PBH) lifetime, for a given initial 
mass. In particular, initial mass of PBHs evaporating to- 
day can be 10^ — 10^° g rather than w 10^^ g as predicted 
by GR. 

For the PBHs having small masses, there are astro- 
physical constraints on their abundance, based, e.g., 
on studies of extragalactic photon and neutrino back- 
grounds. These constraints give, as usual, the informa- 
tion about primordial density perturbations (we assume 
that PBHs form from these perturbations). For an ex- 
traction of this information one must know the evolution 
of these perturbations in radiation era. In the recent 
work by Cardoso et al [l^l it had been shown that the 
density perturbations with short wavelengths are ampli- 
fied during horizon re-entry. The magnitude of this en- 
hancement depends, clearly, on a scale of the density per- 
turbations. The smaller is the scale, the earlier the per- 
turbation crosses horizon, and, if comoving wave number 
k is larger than some critical value fee, this crossing hap- 
pens at high energy regime. The straightforward calcu- 
lation of the enhancement factor, in the region of scales 
which are relevant for PBHs with small masses, evapo- 
rating near today, is quite difficult, even numerically, due 



to a very complicate machinery of cosmological pertur- 
bation theory in braneworld cosmology. 

In the present paper we study the dependence of the 
enhancement factor on the comoving size of the density 
perturbations. We carried out detailed numerical calcu- 
lations of gauge invariant amplitudes of curvature per- 
turbations as functions of the scale factor and the cor- 
responding enhancement factors. We found the approx- 
imate solution (of the equations for perturbation ampli- 
tudes), describing the time evolution of the amplitudes 
near horizon crossing. According to this solution, the 
magnitude of the enhancement factor doesn't depend, in 
the high energy region, on the comoving scale. Using 
this conclusion it is possible to calculate the enhanced 
perturbation amplitudes for arbitrarily small scale. 

The plan of the paper is as follows. In the sec- 
ond section the equations of perturbation theory in RS2 
braneworld cosmology which are necessary for curvature 
perturbation calculations are given. In the third section, 
the approximate solution of these equations in the high 
energy limit is suggested. In the fourth section the main 
relations characterizing the PBH evolution in the RS2 
braneworld are briefly reviewed. The scheme used in the 
numerical calculations is presented in Sec. [V] The re- 
sults of the paper and conclusions are summarized in the 
last section. 



II. SCALAR PERTURBATIONS IN RS2 MODEL 



A. Braneworld cosmology in RS2 model 



More than ten years ago, in works [23|-|26|, exact cos- 
mological solutions in the braneworld had been obtained. 
It was shown also pvl], for the case when the bulk is 
empty, that five-dimensional geometry of all these cos- 
mological solutions is the well-known 28| Schwarzschild- 
AdS (Sch-AdS) spacetime (i.e., the spacetime with 5D 
black hole geometry), having the metric 



-h{r)dT^ 



(1) 



Here, rfS^ is a metric of a unit 3D sphere, plane or hy- 
perboloid (for K = +1,0, —1, respectively). 



M 



(2) 



K is the curvature of the horizon, M is the mass param- 
eter of the black hole at r = 0, ^ is the AdS curvature 
radius. The most natural physical interpretation is that 
a cosmologically evolving brane is moving in this space- 
time, while for an observer on the brane this motion will 
be seen as an expansion of the universe. If the brane tra- 
jectory is given by equations r?, — T{t), ri, — a{t), where 
t is the proper time of the brane, the induced metric on 
the brane becomes 



-df +a^{t)dT. 
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which is the metric of the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) spacetime. 

The parameter M is unknown, but the value of it can 
not be too large, for the braneworld scenario to be con- 
sistent, e.g., with nucleosynthesis data [l^. We suppose 
that M — and shall consider below only this particular 
case. Further, we shall consider the spatially flat brane 
only, i.e., K = 0. Introducing a new spatial coordinate z 
by relation z = £/r, the metric ([T]) with M = 0, A' = 
becomes conformally flat, 

(5)rfs2 ^ _ ^_^^2 ^ ^^2 ^ s..dx'dx^) . (4) 
On the brane, the connection of r and t is given by (T3| : 



n = T{t), T=-Jl + p(- 

ay \ a 



(5) 



which is necessary for static solutions to exist in RS2 
model. The 4D gravitational constant becomes 
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a 



(12) 



At last, the bulk Einstein equations. Gab — —^59ab, 
give the relation between the 5D cosmological constant 
and the AdS curvature radius, A5 = — 6/£^. 

Further, the tensor H^jy in Eq. ([9]) is given by the 
expression 



(13) 

and is the limiting value on the brane of the electric 
part of the bulk Weyl's tensor. The latter term is, in the 
effective Einstein's equations, an external source, with an 
energy-momentum tensor T^f^, defined as 



and a z-coordinatc of the brane is Zf, — £/rh = £/a. 

The main hypothesis of any braneworld model is that 
the string theory predicts Einstein gravity in the bulk, 
i.e., the equation 



Gab — hITab 



(6) 



takes place. In our case, the bulk energy-momentum ten- 
sor has the form 

Tab ^ -^Gab + S'XSBS^^S{y - yb), (7) 



where A5 is the bulk cosmological constant, K5 is a 
5D gravitational coupling constant, S'^^ is an effective 
energy-momentum tensor for the brane, yt is the brane 
position (transverse coordinate of the brane). In the 
Gaussian normal (GN) system, = 0. The tensor S^j^^ 
consists of a brane tension term and the matter energy- 
momentum tensor T^i/, 



(8) 



Using junction condition [29|, one obtains the effective 
4D Einstein equation on a brane (30| : 

(*)G^. = + k'T^. -f 4Ii^,u - f^-- (9) 

In this equation, the quantities k and A4, which are 4D 
gravitational coupling constant and 4D cosmological con- 
stant, respectively, are given by relations: 



A4 = i f A5 + 

2 I ^ 6 
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(10) 



In AdS bulk, A5 < 0. In addition, we will use the RS fine 
tuning condition: 



A. 



5 J2 



(11) 



-F T 



0. 



(14) 



In the case, which we consider in the present paper, this 
tensor is equal to zero because the Weyl tensor, Cabcd, 
vanishes for an AdS bulk. 

The n^^-tensor term in the effective Einstein equations 
(components of this tensor are quadratic in p) leads to the 
following modification of the Friedmann equation (k^ — 
SttG): 



H' = 



8ttG 



(15) 



Deriving this formula, the fine tuning condition, Eq. 
PT|). and equalities K = M — Q in Eq. ^ are used. 
The solution of Eq. ([15]) for a radiation-dominated state 
{p = p/3) on the brane is (see, e.g. 
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a{t) — a 
H{t) = 



1/2 
2t + tc 

3 



32TTGt(t + tc) 



(16a) 
(16b) 
(16c) 



where tc = ^/2. The conservation equation has the same 
form as in 4D case: 



p=-ZH{p + p). 



(17) 



As one can see, at late times (i.e., at low energy den- 
sity), the well-known relations of 4D cosmology are re- 
covered. The time dependence of the horizon mass in RS 
model is (we put c = 1) 



Mh{t) = 



_ Anp _ U'^it + tcf 



3i73 G{2t + tcf 



(18) 



The transition between the so-called high energy (HE) 
and low energy (LE) regimes happens at the "critical" 
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epoch, at which Hi = 1, and horizon mass at this time 
(this happens at t = = t^/ a/2) is 



Mh{tc/V2) w 5 X 10^''^ 



0.1 mm 



(19) 



The critical vahie of the comoving wave number, kc, 
which corresponds to this critical epoch, can then be 
written using the known relations between the horizon 
mass and the comoving wave number (see, e.g., [32] ) as 



M, 



eq 



3 X 10^° Mpc"^ 



-1/12 



9*eq 
1/2 



0.1 mm 



VlOO/ 



,(20) 



where keq, Meq 



and 



are wave number, horizon mass 
and effective number of relativistic degrees of freedom 
corresponding to the time of matter-radiation equiva- 
lence and (7,c is the number of relativistic degrees of free- 
dom corresponding to the critical epoch. 

The main parameter of the model, £, can be con- 
strained by Newton's law tests in table-top experiments. 
The most recent results [s^, [s^] give the following limit, 
which is very important for cosmological implications of 
the model: 



e < (0.015 - 0.044) mm. 



(21) 



The corresponding constraints from astronomical obser- 
vations are somewhat weaker (see, e.g., [3^). 



B. Scalar perturbations 

The case when AI = in Eq. ^ corresponds to a 
pure AdS bulk spacetime. It is known that in this case 
a study of cosmological perturbations in the bulk and 
the brane is greatly simplified. It was shown in 0, [1, [l3| 
that a solution of the perturbed 5D Einstein equations in 
a vacuum AdS bulk, having only metric perturbations. 



^^MGi = 0, 



(22) 



can be reduced to a solution of the evolution equation 
for the "master variable" fl (which depends only on co- 
ordinates in the 2-dimensional orbit space, i.e., on t, z) 
whereas all gauge-invariant metric perturbations in the 
bulk are written in terms of this fl. 

In Poincare coordinate system (used above, in Eq. ([4])) 
the wave equation governing the evolution of the master 
variable in the bulk (the master equation) is 





d^n 3 dfl 




Q^2 


dz^ z dz 


a- 



k^]n = 0. 



(23) 



Here, and everywhere below, we work with Fourier trans- 
forms (with respect to the a;"s) of ft and all the pertur- 
bation functions. 



The important boundary condition for Q. can be ob- 
tained from Israel's junction conditions [l^. These con- 
ditions take the simplest form in a GN coordinate system, 
in which the bulk metric is 



(24) 



The perturbed 5D metric in this system is given, in gen- 
eralized 5D longitudinal gauge, by the expression 



/ -n^{l + 2A) 



9AB 



V 




nA 



V 





(1 + 2n)s, 





1 





-2A 



vv 



(25) 

All quantities in Eq. (j25p and, in particular, n, a, are 
functions of GN coordinates t, y. On the brane one has 
Tih — I, ah — ait, y — 0). The functions a(y, t) and n{y, t) 
are known from the solution of the Einstein equations in 
GN coordinate system [25j. Scalar quantities A, Ay, Ayy, 
TZ in Eq. ((25|) are gauge invariants. The formulas relating 
the derivatives in two coordinate systems are given by 




d_ _l I _^ad_^ 
dy a \ adr 



d_ 

di 



(26a) 



(26b) 



Using the expressions for the junction conditions |36| . 
we neglect in them the terms with anisotropic stress per- 
turbation in the perturbed energy-momentum tensor for 
matter on the brane and, correspondingly, all terms con- 
taining the brane bending scalar £,{t,x^) (describing the 
perturbed position of the brane) in the expression for the 
perturbed extrinsic curvature tensor. In this approxima- 
tion, one can introduce the following notations: 



$ = A, * = -TZb, 



(27) 



having in mind that these gauge invariant perturbations 
of the bulk metric coincide, on the brane, with lapse and 
curvature perturbations in the conventional 4D cosmo- 
logical perturbation theory. 

Junction conditions give the expressions for matter 
perturbations {Sp, Sq, dp) on the brane through the lin- 
ear combinations of gauge invariants and their derivatives 
and, therefore, through the master variable and its 
derivatives. Using these expressions, one can obtain, for 
il, a boundary condition on the brane expressed through 
the g aug e-invariant quantity A (defined below in Eq. 



dy 



1 



1 



6pa 



0. 



(28) 



Considering the perturbed effective Einstein equations. 



(29) 
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one can parameterize the perturbations of Ef^^ in the form 
MM: 



6£^ 



(30a) 
(30b) 



y = e 



K^SpsYS^ + k^SneY^), (30c) 



(31) 

treating the trace free tensor SS/^^i, as an additional fluid 
source term in Eq. (1291) with a radiation-like equation of 
state. This assumption is in full analogy with a case of 
the tensor T)i,y, where one has 



STo" = ~5pY, 

ST," = -kY.Sq, 

ST,^ = SpYS,^ +k^dTrYi^. 



(32a) 
(32b) 
(32c) 



The perturbations of the Weyl fluid, {Sp£,Sq£,6TT£), 
transfer effects of the bulk metric perturbations (effects 
of "KK degrees of freedom" ) to the brane. 

The solution of the perturbed equations is a gener- 
alization of the results of standard 4D cosmological per- 
turbation theory. The corresponding formulas are de- 
rived in (22! (in the approximation 6n = 0). These for- 
mulas express gauge invariants $, in terms of the gauge 
invariant matter perturbation variables A (which is a 
density contrast in the comoving gauge) and V (which 
is a peculiar velocity in the longitudinal gauge), as in the 
4D perturbation theory. These invariants are given by 
the relations (in the longitudinal gauge): 

pA^ Sp-3H5q, a{p + p)V ^ -kSq. (33) 

There are two differences from the 4D theory: the formu- 
las include i) 0{p/a) corrections and ii) KK corrections, 
i.e., the terms proportional to Sp£,5q£ and Stte. These 
latter terms can be expressed through the master vari- 
able [H: 



6p£ 
5q£ 
Stt£ 



k'^n 





tl- 








^3K^a'^ 


a 


1 / 






3f2 








a 



(34a) 
(34b) 



p)n' 



(34c) 



where the prime denotes dy and the dot denotes dt- 

Using the results of [2^ and Eqs. (|34p one can easily 
obtain the ordinary differential equation for the gauge 
invariant A. In the approximation = p/p — 1/3, w — 
p/p = 1/3, one has 



A + HA- 



4p 
a£2 



18p2 



4fc4 



Qb- (35) 



This equation contains the term which is proportional to 
fib, in the right-hand side. Therefore, this equation is 
connected with Eqs. ([23]) and (f28| . 

Another important gauge invariant is the curvature 
perturbation on uniform density slices. It is defined by 
the relation C = ^ — HSp/ p, where ip is the curvature per- 
turbation. The relation between C and A also contains 



c = 



1 3pa2(3p + 2cr) 

4 ^ 4fc2£V2 



3Ha dA fc2 



4fc2 df] Gia^ 



^b. (36) 



III. THE HIGH ENERGY REGIME 

Studying, in the high energy regime of radiation dom- 
inated era (when, in particular, H dy ~ ~dt), the 
dependence of A and Q on time before horizon crossing, 
by power series methods, and taking into account only 
the dominant growing mode, one can obtain (at leading 
order in krj) the result p^ j: 

A"^ « ^(fc??)', nr~3ialk'^{kr,f. (37) 

Here, a* is the scale factor at time of Hubble horizon 
crossing. In the high energy regime one has 



a = a.(3fcr;)i/3. 



(38) 



The connection between a* and the comoving wave num- 
ber is 



a^, = = Qc 



2-1)1/3^ 
' k 



1/3 



(39) 



i{t,/V2) « 1.25 r!^/' ( ^ 



1/2 



10 



-16 



0.1 mm 



1/2 



(40) 



One can rewrite Eqs. ([57]) in the form: 



A- « ± ( ^ 
27 Va* 



9 3 

al 



-A-] tI- (41) 



As one can see from Eqs. (|^T|) and ([M)) . the value of 
A"^* at Hubble horizon crossing is constant and the cor- 
responding value of Q!^^ depends only on k. 

We are interested in a behavior of A and i7b in a rela- 
tively short time interval, from a = a* up to a < 3a*. As 
the results of the numerical calculations show (see Sec. 
IVI|) . just near a ~ 3a* the fif,, A-values reach maximum. 
At later times the oscillations begin, and amplitudes of 
these oscillations are equal, approximately, to the maxi- 
mum magnitudes of i7f,, A reached at the previous period 
of the smooth behavior. 
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Our key assumption is the following: the growth of fit 
amplitude in the interval (a» 3a^,) can be described, in 
the limit of large k, k ^ kc, hy a. function which does not 
depend on the comoving wave number k. Namely, one 
assumes that 



1 



9fc2 



(42) 



in the asymptotical limit of the high energy regime, k ^ 
kc- The function ffi decreases with a growth of a/a* and 
it is assumed that fn{l) — 1. 

According to this assumption, the time evolution of 
fib, starting from the horizon re-entry, is the same for all 
comoving wave numbers fc, and the fc-dependence of fib 
enters only through initial conditions at a = a*. It may 
be justified as follows. A general solution of the master 
equation is given by the expression 38] 



r^b = — I dmS (m) Zo{mz)e 



(43) 



where Zq is the linear combination of Hankel functions 
and S{m) is the arbitrary function, uj = \Jw? + k"^ . The 
variable m has a physical sense of the KK mass. In the 
high energy regime, when Rl ^ 1, the physical sizes 
of perturbations, at horizon re-entry, are smaller than 
a*/fc = <C I. Correspondingly, k <C a^jl. At the 

same time, it is well known that in the high energy regime 
the contribution to fib from the massive KK modes is, 
in general, significant (in contrast with the low energy 
case), i.e., the characteristic values of m contributing to 
the integral for fib can be much larger than a^^jl. So, 
in the high energy regime, characteristic m-values are of 
the same order as fc-values, or even larger, and the fc- 
dependence of fib can be effectively masked (if ujchar = 

The assumption (|^^ is used below, for numerical cal- 
culations of A in the region k ^ kc- The results of these 
calculations show (see Sec. I VI I for details) that in this 
region of fc, the ratios A(a)/A(a,) are the same for dif- 
ferent fc (for a, < a < 3a*). It follows from here that, in 
addition to p2|) . one can assume that 



A 



4 
27 



/a(1) = 1, 



(44) 



in the same interval of a, a* 3a* . 

Substituting now the expressions and for ilb 
and A in the equations used for the numerical calcula- 
tions, one can see that all of them become independent 
on k (in the high energy regime). 

In the high energy regime, when a ^ i^/^, the following 
useful relation holds: 



H^a" ={ — ] e 

a* , 



(45) 



Using this relation and ansatzes p2|) and (|441) . one ob- 
tains, from Eq. ([55)1 for C: 



27 Xa^f I V a* 



o/a ( — 
3 V a* 



6/a(^ 

a. 



-rJ- 



1 



(46) 



n 



54 \a* y \a* 
At a = a*, one has, as it must be, C ~ 1, if the condition 



-fk 



(47) 



holds. At a/a* — 3, i.e., near the maximum, one has 

Cm a 



a „\ 1 ,e " 1 



If /a(3) ~ /o(3) ^ 0.15, one obtains that (max ~ 6. This 
value is close to a value of the enhancement factor (see 
Sec. EH). 

The equation (I55|) for A, after substituting of Eqs. (HH 
B4|) . becomes 



a* \a* J \a* J 



(49) 



Neglecting in Eq. (|49l) the terms with derivatives, in 
accordance with Eq. (|T7)) . one obtains the approximate 
result 



/a(- 

a* 



/.(- 

.a* 



Analogously, from the equation 
condition one obtains 



-fh - 8/a + 8/0 = 0. 



(50) 

for the boundary 
(51) 



For consistency with Eq. ([5(1)) . the function /si(a/a*) 
must obey the inequality 



'fn 



In 

Oi*: \ a* 



(52) 



This condition is consistent with the similar condition 
(|T7)) and with 

In conclusion, we showed in this section, that the as- 
sumptions and (jU) lead to the independence of the 
ratio A(a)/A(a*) on fc in the interval a* 3a*, in the 
high energy regime. Estimates show that for the mode 
with fc = lOfcc the critical epoch corresponds to the mo- 
ment of time when a = 3a*. Therefore, beginning from 
fc « lOfcc, the interval (a* 3a*) is entirely in the high- 
energy regime. Correspondingly, the asymptotical region 
in which A(a)/A(a*) is independent on fc begins from fc's 
which are larger than lOfcc (say, from fc 30fcc). 



7 



10 



15 



15 



10 



o 




Using the expression for Mh (Eq. (IT8l) l and the relation 
15 tc — £/2, one can see that the equality ([55)1 is consistent 

with the inequality (j55p only if i < ic, i-e., in the high 
energy regime. It means that PBHs which form in the 
high energy regime are 5D black holes. 

A rate of a loss of the PBH's mass, due to the 5D- 
evaporation, is proportional to rj^ for the evapora- 
tion in the brane as well as in the bulk. So, one has 
dMBH/dt ^ Mg\j. Resulting lifetime of the black hole, 
tpjiam is proportional to M'^^ rather than as in the 



5 10 15 

log,„ t [s] 

FIG. 1: Black hole mass Mgu versus the moment of time at 
which it evaporates, assuming the case of RS cosmology with 
£ = 0.1 mm. Labels "tjvs", "teq" and "to" show, correspond- 
ingly, the nucleosynthesis, matter-radiation equivalence and 
present epochs. 



IV. CHARACTERISTICS AND EVOLUTION OF 
5D BLACK HOLES 

The formation and evolution of PBHs in RS2 cosmol- 
ogy had been investigated in [3l|, H^-H^] . 

Supposing that the braneworld PBHs localized on the 
brane are represented by the 5D Schwarzschild solution 
[3, [iHl for the metric, one obtains, for the induced 4D 
metric on the brane, the expression 

(53) 

which does not coincide with the 4D Schwarzschild met- 
ric. Correspondingly, the relation between PBH mass 
Mbh and radius. 



ds^ = — 


1 - 


(7)1 


dt^ + 


1 - 


(t)! 
















(54) 



is different from the analogous relation in the 4D case 
(throughout this section, we will use, following (s^, the 
convenient notations, in which M4 is the Planck mass, 
£4 — M^^ is the Planck length, t4 = £4 is the Planck 
time). 

It follows from Eq. that if PBH's radius at its 

formation is smaller than AdS radius £, the following in- 
equality for PBH's mass holds: 



Mbh £_ 
M4 ^ £4 ' 



(55) 



We assume, as usual, that PBHs form with masses equal, 
approximately, to the horizon mass at the time of 
formation, 



evap J 

4D case. 



^evap 

IT 



£i 



e 

M4 



(57) 



In this formula is the time of the onset of evaporation, 
MBH{tc,tevap) IS the PBH mass at t = which cvapo- 




Q 20 -^^ 



-2 

iog,„ m) 




log^„ k[Mpc-^] 



BH 



(56) 



FIG. 2: The dependence of horizon mass on cosmic time 
t (upper panel) and comoving wave number k (lower panel). 
Solid curves show the case of RS cosmology assuming that 
£ — 0.1 mm while dashed curves are for the case of standard 
4D-cosmology. 
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rates at t = t^vap (one assumes that tc ^ tevap)- Here we 
assume, following [l^l, that in the relatively short period 
of time from the PBH formation, t;, up to the end of the 
high energy regime, tc, black hole does not evaporate, but 
increases its mass due to the accretion. The increase of 
mass due to the accretion is determined by the equation 
[13, im dM/dt ~ qM/t {q is the (unknown) parameter of 
an efficiency of the accretion, < g < 1). It is assumed, 
for simplicity, in a derivation of Eq. (j57p that at t = tc 
the accretion process ends completely, giving place for 
the pure evaporation. 

One can check that if AdS radius is too small, the con- 
dition for 5-dimensionality of PBHs, rg < £, can not be 
satisfied. Comparing the mass-lifetime relation ((57)) with 
the expression for the radius (j54l) . one can determine, for 
a given lifetime, the minimal possible value of £, given by 
the relation 



evap 



1/3 



(58) 



If, e.g., PBH evaporates today, tevap ^ to ^ 10^'' s, one 
has = lO^o^. 

The BH mass at the onset of the evaporation (< = tc), 
if the age of the black hole is equal to the age of the 
Universe, is 



MBH{tc,to) = Af*^(to) « 5 X 10^ 



-1/2 



£ > W£i. (59) 



and, in general, for PBH evaporating at time t 

£ 



M-^„{t) « 5 X 10^ 



-1/2 / ^ \ 1/2 



> 



1/3 



(60) 



The dependence of A/^^ on i for ^ = 0.1 mm is shown 
in Fig. n 

One should note, in a conclusion of this section, that 
if the accretion efficiency is not small, the initial masses 
of the PBHs are smaller than its masses at the onset 
of the evaporation. And, even without any accretion, 
initial masses of the 5D PBHs are, for the same values 
of total lifetime, much smaller than initial masses in the 
standard cosmology. It means that in the 5D case the 
known astrophysical constraints on the PBH abundance 
correspond to primordial perturbations on smaller scales. 

The dependence of horizon mass on time and k in brane 
cosmology is shown in Fig. [21 It is seen that PBHs with 
mass less than ~ 10^^ g are produced in the high energy 
regime. The corresponding comoving wave numbers are 
larger than 10^^ Mpc~^. If £ « 0.1 mm then PBHs evap- 
orating today correspond to fc > 10^^ Mpc~^ ~ lO^kc- 



t ^ 



At 



.biane 



II 



tt 



6 = 1 




FIG. 3: (upper panel) The computational grid on which we 
solve the system of equations psp and (|64p with boundary 
condition (|28p . The physical brane is at ^ = ^jv = — 1 while 
the regulatory one is at ^ = = 1- (lower panel) An example 
of computational result for master variable f2 (in this case, we 
took k — 30fcc and A*' — 64, so brane is at ^ = ^64 = The 
grid over n is shown to be homogenous, however, note, that 
the actual ^n-grid is inhomogeneous due to Eq. 



NUMERICAL SCHEME 



It follows from Sec. IIVI that in brane cosmology PBHs 
having, at formation, relatively small masses (< (10^ — 
lO^O) g) and, in particular, evaporating near today, had 
been produced long before the critical epoch, tjorm ^ tc- 
The corresponding comoving sizes of perturbed regions 
are also small, fc^^ - (10^^^ - lO^i'^) Mpc. Therefore, 
it is practically important to determine the enhancement 
factors for rather large values of fc, fc > (10^ — 10'')fcc. 
The straightforward numerical calculation of these fac- 
tors for such large k are quite difficult and unreliable, 
but, luckily, the numerical calculations for moderately 
large k, k ^ (10 — 30)fcc, show the flattening and the pos- 
sible saturation of the dependence of the enhancement 
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factor Q on k. We argue now that, really, Q{k) does not 
depend on k in the limit of very large k, k ^ kc- 

For the numerical solution of the system of equations 
(1231) and ([55)1 with boundary condition (1^51) . a pseudo- 
spectral calculation method was employed. Such meth- 
ods are often used in the tasks of hydrodynamics and a 
detailed description can be found, e.g., in 

To be able to perform a spectral transformation over 
the set of Chebyshev polynomials we do a following 
change in the variables: 



2z- 



(Zreg + Zb{t)) 



Zb{t) 



i<e< 1, 



(61) 



(62) 



where t has the meaning of cosmic time on the brane and 
is related to t by 



dr 
'dt 



a{t) 



(63) 



while Zreg IS the position of the regulatory boundary (ar- 
tificial cutoff that is introduced to make a computational 
domain finite, see, e.g., [H, E^l). 

The equation for f2 (1231) is rewritten in new variables 

as 



9*2 

where 



dm. 



K, 



'dt 



drKdr 



im = 0, 

(64) 
(65a) 



Zb{t) 



dr^ [dTj 



= 2 



6 



,9^ 

z{t,^){Zreg - Zb{t)) 

d^l ■ 



(65b) 
(65c) 

-2 

(65d) 
(65e) 



Further, a new variable x related to time derivative of 
57 is introduced [i^ to reduce the task to two first-order 
equations: 



= X-Kt^^ = F{x,n'^:. t,0, 



dt 



(66) 



9x 

dt 



-K: 



dHl fdKt^ 
W ^ V dt 



Kt X-K, 



dn 

9? 



dn 

-Kn = 



To solve this system using difference method, the 
transformation of all quantities over set of Chebyshev 
polynomials is done for the ^ (and x) ~ axis. This is 
done at each time step so that system of partial differen- 
tial equations (l66l |67| reduces to the system of ordinary 
differential equations. 

Thus, at each point (ip,f„), the following quantities 
are known: 

x,n,n'^,n'l^,F,G, 

and also known are Chebyshev transforms 

Xni nn, {n^)n, (f7^^)„, Fn, Gn- 

The grid (see Fig. |3] for illustration) based on Gauss- 
Lobatto points is used here, 



/TTn 

= COS ( — 



n = 0,l,...,iV, 



(68) 



because it allows to perform fast Fourier transforms 
(FFTs) between the set of values of any variable (e.g., 
n) in Gauss-Lobatto points and its Chebyshev compo- 
nents r2„. Chebyshev transforms of derivatives, such as 
(ri^)„ and (n'^^)n, are also easily obtained using recur- 
rence relations from the Chebyshev components of the 
function (see [11] for details). 

Equations that are actually solved on each time step 
are: 

^-Fn{t); ^ = G„(0. (69) 

For points on the brane, A(ip) is also evaluated at each 
step and Eq. (I35p is solved using a finite difference 
method (in our case, a 4-th order Adams-Bashforth- 
Moulton scheme). 

The boundary conditions are imposed on the values of 
the highest two components of the master variable, 17 at 
and r^Tv-i. This is done by demanding the following: 



N 



f7^(-l)-^aX(-l), 

n=0 
N 

^(i) = ^fi„T;(i), 



(70a) 
(70b) 



71=0 



where T'„(^) is the n-th order Chebyshev polynomial. 
The value of 57^ (1) is assumed to be zero (condition on the 
regulatory brane) while the value on the physical brane, 
r2^(— 1), is related to other quantities by the boundary 
condition ([^5]) which can be expanded as 



dn\ 



a(zr^g-zt{t)) a dr ^ yi+TfTp" ^ 

_ (71) 
Following the approach of [4a], we also use the following 
(67) additional condition: xn = Xn-i = 0. 
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gree of enhancement of the perturbation amplitudes: 



A, 



eff 



^GR 



^5D 



Q5D = 



A 



5D 



^GR 



= QeffQe 



(72) 

In a case of the effective theory (fif, = 0), the en- 
hancement reaches an asymptotic value, Qeff ~ 3, at 
k ~ lOOfcc. However, the direct calculation of Q^d (or, 
equivalently, Qg) for very large wave numbers A: 3> fc^ is 
not easy, due to a quite complicate behavior of fl in the 
bulk (see Fig. [3] for an illustration: the larger value of k, 
the more frequent are the oscillations in the bulk). Due 
to limitations of computing resources, we have been able 



a -2 



o 




-0.2 0.0 



0.2 0.4 0.6 0.8 



FIG. 4: The result of the numerical calculation of master vari- 
able n for k = 3fcc. Upper panel: the illustration of 
(here, A'^ — 32, units over time and f2-axis are arbitrary); 
lower panel: the value of Q on the brane, ili,. The normaliza- 
tion is given by = 1 for a <g a,. 



VI. RESULTS AND CONCLUSIONS 



The main results of the numerical calculations are 
shown in Figs. |4]- [H Figs. SI [5] show, for a given value 
of k, k = 3kc, the evolution, near horizon crossing, of 
three variables: fit (Fig. 01), A (Fig. O upper panel) and 
C (Fig- 13 lower panel). It is clearly seen that all three 
variables rise with a, up to a ~ 3a*. It is also seen from 
Fig. [5] that the corresponding rise of A and C in GR is 
weaker, and, as a result, we have an enhancement. It 
follows also, from Fig. [Sj that the enhancement is not 
zero in the approximation flj, — 0, when there are no 
KK corrections in the equations of the 5D perturbation 
theory. 

In Figs. [6l [7] it is shown how the evolution curves for 
and Qb change with an increase of k. It is seen, from 
Fig. [51 that the enhancement grows with fc, but there is a 
clear tendency of a slowdown of this growth at k > lOfcc- 

Following [2^ we define the factors that show the de- 



0.5 



0.0 



o -0-5 



-1.5 




-2.0 

-0.2 0.0 0.2 0.4 0.6 0.8 1.0 
log,„ (a/a,) 

-0.2 0.0 0.2 0.4 0.6 0.8 1.0 



O 




0.2 0.4 0.6 

login 



FIG. 5: The calculation of the density contrast and curva- 
ture perturbation on the brane for k — 3kc for the cases of full 
5-dimensional calculation (curves labeled "5D"), effective ap- 
proach (approximation of Q,i, — 0, curves labeled "Eff") and 
General Relativity, i.e. standard cosmology ("GR"). Note 
that the value of a* is the same for "5D" and "Eff" cases, 
but is different for "GR" because of the different evolution 
of the background quantities. Upper panel: Comoving den- 
sity contrast A as a function of the scale factor, normalized 
to ^ = 1 in super-horizon regime. Lower panel: Curvature 
perturbation ^ calculated using the same three approaches. 
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O 




0.2 0.4 0.6 

logio (a /a*) 

FIG. 6: The result of the numerical calculation of for dif- 
ferent values of k. 



o 




-16.5 



-16.0 -15.5 

log a 



-15.0 



FIG. 7: The result of the numerical calculation of fit for 
different values of k, normalized to ("fe = 1 in the super-horizon 
regime. Arrows show the horizon crossing time (fc — aH) for 
each mode. 



to make direct calculations in 5D case only for a limited 
range oi k < SOfcc- 

In order to study the PBH production for masses 
Mgjj{to) ~ 10^ g (such PBHs, as we have seen in Sec. 
rvl evaporate near today, if the value of £ is close to its 
upper bound (j2ip ). we need information about cosmolog- 
ical perturbations for k > lO^k^. To perform calculations 
for such large wave numbers, we have used an approxi- 
mate approach according to which J7f,(a/a*) has the same 
form (and is given by Eq. (j42p ) for all large wave num- 
bers (see the discussion in Sec. IIII|) . In these calculations 
we have used, for the required function /o(a/a*), the cor- 



responding function obtained from the direct numerical 
calculation for k = 30/cc- Using this approach, we have 
calculated the enhancement factors for large values of k 
{k > soke). The results of the calculation are shown in 
Fig. 1 

In summary, we stressed in this paper that, in RS2 
brane cosmology, the PBHs of relatively small mass (the 
concentration of which in space can be constrained by 
cosmological arguments) form in the high energy regime, 
and the corresponding comoving wave numbers are very 
large, k ~ (10^ — 10^)kc- We thoroughly studied, by 
numerical methods, the evolution of scalar perturbation 
amplitudes (those needed for calculations of the PBH 
production) near horizon crossing, for a wide range of 
comoving scales. We confirmed the main conclusion of 
[2^ according to which amplitudes of the curvature per- 
turbation get enhanced after horizon re-entry (before a 
beginning of the oscillation phase). We developed an ap- 
proximate phenomenological approach for calculations of 
the perturbation amplitudes for very small scales, where 
the direct numerical methods are powerless. We argued, 
using this approach, that in the asymptotic limit of high 
energies, the enhancement factor is constant as a func- 
tion of the perturbation scale. We presented details of the 
numerical scheme (based on the pseudo-spectral method) 
which is used for a treating of scalar cosmological pertur- 
bations on the brane and in the bulk. 




5D/GR 



Eff/GR 



5D/Eff 



12 3 4 

logio(kAc) 



FIG. 8: Enhancement factors that show the degree of increas- 
ing of the perturbation amplitude after horizon entry. From 
bottom to top, curves show the enhancement of the ampli- 
tude of 5-dimensional calculation compared to the effective 
one, effective theory compared to General Relativity result 
and 5-dimensional calculation compared to General Relativ- 
ity result. 
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